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ABSTRACT 

Using the linear sigma model to describe quark-pion interactions, we 
compute polarization asymmetries in quark fragmentation. We show that 
the effects of transverse quark polarizations appear in the correlation be- 
tween the two leading pions in a jet produced by the fragmentation of a 
quark. Such asymmetries provide a window to the nature of chiral symme- 
try breaking in QCD. 



1. Introduction 

Measuring tiie polarization of final state particles has been of continued interest in high 
energy particle physics. Such observations facilitate the measurement of the form factors 
which govern particle interactions and provide a probe of the symmetry breaking mech- 
anisms present in nature.affl One interesting idea is to probe the polarization carried 
by a quark or gluon of QCD through the distribution of its fragmentation products in a 
final state jet, first proposed by Nachtmann.aa In particular, Nachtmanna showed that a 
three particle correlation in a jet can probe the helicity of the parton initiating the jet. 
This idea was rediscovered by Efremov, Mankiewicz and TornqvistS, and they named the 
correlation the 'handedness' of a jet. 

In collaboration with HeppelmannQ, we showed that, in a similar fashion, a two- 
particle correlation can be used to probe the transverse polarization of the quark initiating 
a jet. This is a novel idea, and there is at present no known experimental information 
on this two-particle correlation. In this paper, we will perform some very simple model 
calculations of the transverse spin dependence of fragmentation. We will use a linear sigma 
model of pions coupled to quarks, somewhat in the spirit of Georgi and Manohar.H We will 
assume that this model, taken in the lowest relevant order of perturbation theory, is an 
approximation to the long distance dynamics of QCD. Within the model, we will calculate 
the spin-dependent part of the fragmentation of quarks to a two-pion state, q ^ rnv + X . 
The spin dependence shows up as a correlation between the plane of the two-pion system 
and the transverse spin vector of the quark. 

The importance of our calculation is that it shows that there is consistency between 
the symmetries of QCD and a nonzero transverse-spin-dependence for fragmentation at the 
leading- twist level. Indeed, the large effect we calculate shows that there is no suppression 
of the spin-dependence. The spontaneous breaking of chiral symmetry (and hence a non- 
zero pion mass) is essential to our calculation. Moreover some kind of nontrivial phase 
or interference is also essential. Since we are in a strong-coupling regime, this does not 
preclude a large effect. Indeed, the size of the effect is very large — 50% or more. In our 
simple model, the interference is between two-pion production in the continuum and at 
the a resonance. 



After reviewing the sigma model and the methods for computing decay functions, 
unpolarized quark fragmentation to a a is calculated. This is used to normalize the polar- 
ization asymmetry. We then compute the polarized decay function to a two-pion state, and 
establish a nonvanishing asymmetry. This then demonstrates how the correlation between 
the two pions reflects the spin of the quark. 

Finally we present some numerical calculations. 



2. Sigma Model, Feynman Rules and Decay Functions 

2.1. Sigma Model Lagrangian 

We modify the sigma model lagrangianll^ to use quarks instead of nucleons: 
1 
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where 



(d^a)' + {8^3)' +iqjt'a^q-gq{a + ir-7i^i)q-V{a' + g') (2.1) 



V = 4 f"^' + *') + ! ("' + *')' 



The Lagrangian possesses an SUji{2) x SUi{2) symmetry, and the term linear in the a field 
that is normally used to provide explicit breaking of the chiral symmetry is not relevant 
for our purposes. 

The sigma field is an isosinglet while the pion field is an isovector whose components 
may be written 

(2.2) 




By using the isospin invariance of the pion field we can define the charged and neutral pion 
fields, 



TT 



±- 1 .^ ^.^^ ^0 



' V2 
The quark field is an iso doublet with three colors 



TTl =F ^7r2) , Tl = TTS. (2.3) 



d ) ■ (2.4) 



The {T°')ij are the 2x2 Pauh matrices, 



.i_/oi\ 2_/o-n 3 n 



We use an index a = 1, 2, 3 for the pion fields and indices z, j = 1, 2 for the quark fields. 

Because of the wrong sign mass term in Eq. (2.1), there is spontaneous symmetry 
breaking. We choose the a field to have the nonzero vacuum expectation value: 



(0|(7|0) = ^, (0|7r|0) = 0, v^^fi^. (2.6) 

After the transformation a ^ a + v, the quark interaction becomes 

-gq {a + ir- 11^5) q — > -gvqq - ^^ (a + if ■ ^75) q, (2.7) 

and the remaining n and a interactions become 

_y (^2 ^2^ ^ _ 2 2 _ ^3 _ A^4 _ A .4 _ A^2^2 _ ;^^^^2_ ^3.8) 

K ' J ^^ 4 4 2 

Note that the pion field is massless, the a has a mass of v2a* ^^^ the quarks have gained a 
mass of gv. If we take the up quark mass to be around 300 MeV and the vacuum expectation 
value of the sigma field for the spontaneous symmetry breaking to be f = /tt = 92 MeV, 
then ^ ~ 3. Typically, the sigma mass is taken to be aroundEJJ 600 MeV. The couplings g 
and A are large, as is appropriate for hadronic interactions in the non-perturbative regime. 

2.2. Feynman Rules and Decay Functions 

Among the Feynman rules obtained from the Sigma Model lagrangian after spontaneous 
symmetry breaking are the vertex factors, as displayed in [big. 2.1|. To perform a calculation 



involving charged pions it is useful to introduce isospin "polarization" vectors: 

£« =(l,±z,0)/v^, 4 = (0,0,1). (2.9) 

An incoming state gets an £*and an outgoing state gets an e*. A color average over the 
initial state is used. 

In our model, we will compute an approximation to the spin-dependent fragmentation 
of a quark with helicity h and space-like transverse spin vector Sy. The total spin will 
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Figure 2.1: The relevant Feynman rules obtained from the Sigma Model 
lagrangian. 

obey the constraint < /i^ — Sy < 1. (Note that with the usual metric, Sy < 0.) In 
the definition of the fragmentation functionaZIo there is a trace with the following Dirac 
matrices: 

G = ^ (7+ + /17+75 + 7+75 ^r) • (2.10) 

The light-cone component, manifest in the overall factor of 7^, provides the projection 
on the leading twist part of the fragmentation. Eq. (2.10) gives the Feynman rule for 
the cut eikonal vertex that defines the fragmentation function. The 7^ term gives the 
unpolarized fragmentation. The helicity dependent term h'j'^'y^ will not contribute since 
we only measure two particles in the final state; a helicity dependent asymmetry only 
occurs with the three-particle correlation called "handedness" .Q'&Q 



Figure 2.2: The fragmentation of a i/' into hadron A. 

The Feynman rules for fragmentation functions were written down in [Ret. 12| , with 
the polarization dependence given in [Ref. 7|Ref. 8| . The distribution of a single measured 
hadron A of momentum p in the fragmentation of a quark is a function of the longitudinal 
momentum fraction z of the hadron and of the transverse momentum ky of the quark 
relative to the hadron. It has the formcl 



Da/^ (2, kr) 



dk- 



:^A/^{k^,k-p) 



(2A1) 



where k'^ = p^ j z, and ^ aj^I)^ •, k ■ p) is the inclusive two-point function, 
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^A/i.{k\p-k)=- / d'xe''-^TvG{0\^{x)al^,^{p)aontA{p)i^{0)\0) 



(2.12) 



with the Dirac matrix G being defined in Eq. (2.10). The overall factor 1/6 represents an 
average over quark color and spin. This is represented diagrammatically in |Fig. 2.2 . The 
factor CiQ^iAiP) ^ont a{p) is the number operator which counts how many hadrons of type 
A and momentum p are in a state. Note that the transverse momentum of the hadron 
relative to the quark is — zky. 

The usual decay function is obtained by integrating over transverse momentum: 

(2.13) 



DA/i, (z) = d \<^tDa/^ {z, kr) . 
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It has no dependence on the quark spin (/i, s). To interpret the fragmentation function as 
a probabihty would require an extra factor, but since we wiU be computing polarization 
asymmetries, this factor will cancel between the numerator and the denominator, so we 
will ignore it. For further details the reader is referred to [Ret. 12| and [Ret. (j| . 



Exactly similar formulae apply to fragmentation to 2 measured particles. We sim- 
ply replace the number operator %-^^j^{p)aoutAip) by the corresponding operator for two 
particles: 

4ut.li (Pl) ^lntA2 (^2) aout^2 (^2) flout^i (Pl) , 

SO that the fragmentation function to two particles is defined by 

'^ ^J 12{2nf (2.14) 

X Tr 7+(0|^a (0, y', yj_) \AiA2X){AiA2X\tl;a (0) |0). 

The polarization-dependent part iaOo 

' ^J 12(27r)'^ (2.15) 

X Tr7+757i ■ Si(0|V'a (0,y",y±) \AiA2X){AiA2X\i^a{Q) |0). 

Ultimately we will integrate over /cj^, the transverse momentum of the quark relative 
to the measured hadrons, and in particular this will imply an azimuthal average over /cj^. 
In this situation, to get transverse-spin dependence will require that there be an imaginary 
part in the amplitudes making up Eq. (2.15) (since the factor of 75 will give a factor i). 
In lowest order for the fragmentation q* -^ nnX, the imaginary part will come from the a 
propagator, as in [Fig. 2^ (c) and (d). 



First, we will compute the fragmentation for q* —>■ aX , followed by the decay a -^ tttt. 
This will give us the total rate for q* -^ nnX at the resonance. Then we will compute the 
spin dependence from interference graphs for q* -^ nnq, with a correct treatment of the 
imaginary part of the denominator of the a propagator 

The momentum of the fragmenting quark is denoted by k and, depending upon the 
situation, the momentum of either the a or tttt pair is represented by p. The Lorentz frame 
chosen is that where the sigma particle (or rather the combined momentum of the two pions 
of interest) has zero transverse momentum. The longitudinal momentum fraction of the 
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Figure 2.3: The four cut diagrams for tlie u* -^ tt+tt u fragmentation 

two-pion system (or a) is given by 2 = p^ /k^ . Generally, these two-pion system variables 
will be written in terms of individual pion variables, 



Z = Zi + Z2, P=Pl+P2, 



^2.16) 



where the individual pion variables will be defined according to the process considered. 



2.3. Decay Width of Sigma 

The imaginary part that we need for the spin-dependence comes from the imaginary part 
of the denominator of the a propagator. At lowest order, this will be given by the width 
of the a. Recent work has appeared where the random phase approximationlld and one- 
loop computational^ have been applied to the linear sigma model to compute the width 
r(o" -^ nn). However, since we are simply demonstrating the existence of a spin observable 
in this paper, we shall not enter the quantitative detail of these references. 



With the cr— TT— TT couphng described by |Fig. 2.1| , the decay width may be expressed as 



r(a^7r7r) = ^^, (2.17) 

where two thirds of the width is due to the decay to charged pions and one third of the width 
is due to the decay to the neutral pion pair. Using m^r = 600 MeV, niq = 300 MeV, and 
f = 92 MeV, we find that r(a -^ nn) = 761 MeV. The relations between the parameters 
which exist due to the spontaneous symmetry breaking, allow us to write the width in 
other forms: 

r{a^..) = ^ = ^. (2.18) 



3. Fragmentation q* -^ aX 

Our goal is to compute the tttt asymmetry near the a pole, so we next find the total 
production rate for fragmentation to a. This rate will provide the normalization for the 
polarization asymmetries. The distribution to be evaluated is represented by the diagram 
in [Fig. 2.3| (b), but with the pion loop removed. The momenta are defined such that 
k represents the incoming virtual quark momenta while p3 and p respectively represent 
the outgoing quark and sigma momentum. (The momenta pi,p2 will represent the pion 
momenta in q* —>■ qnn.) The + momentum fraction carried away by the sigma is given by 
z. Summing over final state quark spins and averaging over initial state spins and colors 
yields the relation 

J (271) 4(fc2-m2) ,^3j) 

5 {k+ - p+/z) S ({k - pf - mfj , 
where the computation is being done in d = 4 spacetime dimensions. There will be no 
spin dependence for a single particle fragmentation after we integrate over the transverse 
momentum (or at least its azimuth). 

The definition of the fragmentation function is that we use a frame in which the a has 
zero transverse momentum, pr = 0. A simple calculation yields 

q2 2 m^J2-zf + z^k^ 

D^/ (z, kr) = ^ ^- ^. (3.2) 

^' 167r3 [m2(l-^)+m2^2 + ^2^2^]2 
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Normally, fragmentation functions are defined with an integral over /cy, the transverse 
momentum of the quark relative to the measured hadrons. These integrals have divergences 
at large /cy which are removed by renormalization (or some other means, like a cut off). 
The dependence on the cut-off or on the renormalization scale gives the usual Altarelli- 
Parisi evolution. For our purposes, however, it will be convenient not to bother with the 
integral. In any event, our model is only appropriate integrated over low kx, where QCD 
is non-perturbative: The fragmentation functions would then be suitable initial data for 
normal perturbative QCD evolution. 

Therefore we will not present the results integrated over /cy. 

To obtain the resulting fragmentation function to tt+tt^, we must insert the decay 
o" -^ 7r+7r^ into the quark fragmentation derivation. This is given by the sigma propagators 
and the pion loop in [Fig. 2.3| (b), which result in Eq. (3.1) being multiplied by a factor 



i{-2i\v) 



(3.3) 



p2 _ j^2 _|_ IfYi^Y 

We are aware that the width of the sigma is in the neighborhood of its mass (using 
rricj = 600 MeV, niq = 300 MeV, and v = 92MeV), but since our main interest is to 
demonstrate a nonzero asymmetry, and its order of magnitude, we just use the narrow 
width approximation, 

2 



P 



2 — m'i + inifjT 



-5 [p^ - ml) . 



(3.4) 



For the total rate, and in the narrow width approximation, the graphs in |Fig. 2.3| other 



than (b) are unimportant, since they are less singular at p^ = m^. Hence the lowest order 



approximation for unpolarized fragmentation to tt+tt pairs is 

327r2 



D 



TT + TT /q 



-D^/q{zi+Z2,'kT). 



(3.5) 



The factor 327r^/3 can be considered as IGtt^ times the branching ratio, 2/3, to 



n ' n 
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4. Polarized Partons and Two— Particle Fragmentation 

We now compute the dependence on transverse spin of the fragmentation q -^ tt+tt^X, 
using the definition Eq. (2.15). Since our asymmetry will be at the a pole, we define the 
sigma momentum as the combined momentum of the two final state pions. 

We consider the decay of an up quark which produces a charged pion pair. At lowest 
order, there are four Feynman diagrams for the u* -^ n^n^u fragmentation, and they are 
shown in [Fig. 2.3| . Topologically, this decay is the same as the decay for the down quark, 
d* -^ n^n^d. Likewise, nonzero asymmetries are also expected from the fragmentation to 
TT^— tt'^ pairs. The quark decays to tt'^tt'^ pairs will produce a smaller asymmetry, because 
the azimuthal dependence is sin (p, and Bose symmetry of the two pions therefore requires 
a zero in the spin asymmetry at zi = Z2- 

4.1. Getting the Polarization Dependence 

All the spin dependences for the lowest order u* -^ tt+tt^w are in the imaginary parts of 
the traces from the cut diagrams. Since the decay functions must be real, combining the 
results from all of the relevant tree diagrams will result in the cancellation of all of the 
spin contributions for the total decay functions unless there is another factor present with 
an imaginary part. Such is the case when we are at a resonance. 

Of the diagrams in [b'ig. 2.3| for u* —>■ nnu, the only ones which can contribute a 
transverse spin dependence to the decay functions are (c) and (d). They have interference 
between the continuum and a a propagator. We next note that these diagrams each have 
one factor of the a-propagator, for which we can write 

-2 T-— F = P^^^ + ^^ (P^ - "^') ' (4-1) 

in the narrow width approximation. The sum of the lowest order decay function amplitudes 
represented by |Fig. 2.3| (c) and (d) can then be expressed as 

a + ib 



23? 



(4.2) 



p2 _ j^2 _ ifYi^Y 

Here a and h are smooth functions of the kinematic variables in the vicinity of the a 
pole. Moreover, as we will see, the spin dependence resides only in the h coefficient. 
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The principal value in Eq. (4.1) implies that the spin-independent term gives a small 
contribution when integrated over a neighborhood of the a pole, so that Eq. (3.5) gives 
the unpolarized fragmentation. Nevertheless, it is Eq. (4.2) that gives the singular part of 
the spin dependence: 

2n5 {{pi + P2f - ml^- (4.3) 

The decay function can then be expressed as the sum of an unpolarized contribution 
and a spin asymmetric piece 



Dn*- 



.nnq (^1, ^2, ^T, ^t) = D^.^^^^ (^1, ^2, ^t) + ADj*^^^^ (zi, ^2, ^T, kr) • (4.4) 

In the narrow width approximation, both the unpolarized and the polarized term have a 
factor 5{{pi +^2)^ — ^a); so that the ratio, which is the spin asymmetry, is well-defined: 

,pol 

I ,61 . ,67. O'l'. S\.'r I 

(4.5) 



AL>jri,^^g(2i,22,sr,kr) 



3^Z),/,(zi + Z2,kr) 
Here, we used Eq. (3.5) for the denominator. 



4.2. The Amplitude Computations 



We now compute the necessary terms from the diagrams in [Fig. 2.3| for the fragmentation 
of an up quark into two charged pions. We denote tt"*" momentum by p2 and the tt" 
momentum by pi, and we use a label {a,h,c,d) to indicate the particular diagram in 
IFig. 2.31 



-2'K~^\vCc ^ ^^^^ + m,) (^ + m,) G (^ + m,) (^- ^2 - m,)] 



[k'^-mf) {k-p2) 



mi, 



{Pl + P2) -ml + iniaT 



(4.6) 



5 ({k-pi -p2f -ml 

The Cc factor carries the isospin factors and for u* -^ n^n^u, Cc = Cd = 2. Graph (a) has 
no singularity at the a pole, so that we do not need to calculate it for our approximation. 
We have already obtained the contribution of graph (b), in Eq. (3.5). Notice that F^ is 
the complex conjugate of Fc, so that we do not need to calculate it explicitly. 
Therefore, we obtain the spin dependent portion of the decay function from 



AD{zi,Z2,ST,kT) 



dk_ 

(2^ 



(Fc + Fd) 



= 3? 



-I (P1+P2) =-^0 



dk 

(2^ 



(4.7) 
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If we equate Fc + F^ from Eq. (4.6) with Eq. (4.2), we find the coefficient b for u* -^ n~^n u 

takes the value 

^_ -2Ti\vfCc5{pi-m'';) 

2 (pi ■ Pi) [A;2 - mf\ 



-] £(pi,P3,s,p) -/img£(pi,p3,p, +) + (s-/c)£(pi,p3,/c,+) 

"2 [P^ + 2p ■ Ps) £ (Pi, P3, s, +) k 

(4.8) 
where ps = /c — pi — P2 is the momentum of the on-sheU final state quark, and p = pi +P2- 
We define e{p,q,r,s) = ^k.x^ivP'^Q^t^^s^ i with the Levi-Civita symbol ei^\^^ obeying the 
following relations in light cone coordinates (LCC), 

gM^^iT = ;p ^M^^os (using 4= normahzation per LCC index) (4.9) 

e {p, k, s, 0) = — f p X A; j ■ s, 

where the convention is set by £0123 = +1- 

Since we will use the fragmentation function after integrating over ky, we will average 
over the azimuth of ky. After this average, the helicity term vanishes because the for- 
mulas no longer maintain a sufficient number of Lorentz vectors to keep the contraction 
with the Levi-Civita tensor asymmetric in its indices. We also apply the narrow width 
approximation provided by Eq. (4.1). We obtain 

An A^ , ^2 2^ ^^Ai>z-(sr X pir) 4 (1 - z^) 

^ ^ 27r2 2pl^ [{kl + m2) zl + ml{l- z.f 

' (A;| + mfj {zfz2ml - 2zizlplj) + (1 - z^) U (1 - z„) ml - Aziz^mfi pj^ 
m^Z2 



(1 - Za) pIt + z\ {k\ + m2)J - 4^2 (1 _ z„yp\^k\ 

(4.10) 
where z^ = z\ + Z2 and z = (0, 0, 1) in rectangular coordinates. 

The main demonstration of this paper is complete. Using the linear sigma model, 
we have achieved a nonvanishing asymmetry in the transverse polarization of a quark at 
leading twist. The nonvanishing component in Eq. (4.10) was maintained through the 
interference of the continuum production of tttt pairs with the sigma resonance. Further- 
more, it is apparent that the existence of Eq. (4.10) is dependent upon a nonzero value 
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for the expectation value of the sigma field, i.e., it is the broken chiral symmetry which 
permits this asymmetry to exist. Consequently, this asymmetry can (theoretically) probe 
the chiral nature of QCD. In the limit that the bare quark mass is zero, nonvanishing 
polarization effects remain and without singularity as long as z is not at its endpoints. 

The asymmetry at the sigma resonance may be estimated by comparing the spin de- 
pendence of (f -^ qim given by Eq. (4.10) with the unpolarized fragmentation given by 
Eq. (3.5). One immediate consequence is that the the correlation of the pions' direction 
with the transverse quark spin is given by the cross product (sy x piy). This is in accor- 
dance with the general theory .0 So, as the pion momentum vector rotates about the quark 
(jet) axis, the asymmetry exhibits a sinusoidal rise and fall as it moves with respect to the 
transverse spin direction of the quark. 

We have used the perturbative approximations, low order graphs and the narrow width 
approximation (F ^ rrio)-, but these are not good approximations. Moreover, the use of 
both quark and hadron degrees of freedom as we have done is a crude model. So the follow- 
ing quantitative calculations must be only considered very rough estimates. Nevertheless, 
they do indicate that the spin-dependence is as large as it can be. 

In [Table 4.1| , we list some values of the asymmetry: the ratio of Eq. (4.10) to Eq. (3.5). 
We have replaced the triple product (sy x pir) by |pir|, so that the numbers represent 
the amplitude of the sinusoidal dependence of the pion production on the azimuthal angle 
about the jet axis. 

Notice that the numbers are large. One might worry that the asymmetries go above 
100%, and that within a lowest order perturbative calculation (from the graphs of |Fig. 27^ ) 



this would correspond to a negative cross section. However, we are working in the neigh- 
borhood of the pole of the a. Therefore, it was essential to use a dressed propagator for 
the a lines, to get a suitable imaginary part and resonance width. This means that our 
model is not totally self-consistent: We have performed a selective resummation of graphs. 
We have further made a narrow width approximation on the a pole: This is obviously far 
from perfect for physical values of the couplings of our model. 

Our approximations are valid in the weak coupling limit, ^, A ^ 1, whence the asym- 
metry is of order ^vA, which is then much less than one. The true strong interactions are, 
of course, strong, and so obtaining an asymmetry above unity from our calculation is not 
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Table 4.1: The asymmetry of Eq. (4.5) is tabulated for various values of 
/cr, zi and Z2 using st xpiT -^ \st\ \pit\- We use the values ma- — 600 MeV, 
92 MeV. 



ruq = 300 MeV, v 



THE TRANSVERSE SPIN ASYMMETRY 


zi = 0.4, Z2 = 0.3 


Zakr (MeV) 


Asymmetry 





1.31 


200 


1.21 


400 


0.89 


600 


0.58 


1200 


0.20 


zi = 0.2, Z2 = 0.15 


zakr (MeV) 


Asymmetry 





0.31 


200 


0.31 


400 


0.34 


600 


0.39 


1200 


0.54 



impossible. All it indicates is that the effect we are calculating suffers from no particular 
suppression. Hence we can expect substantial analyzing power for quark transverse spin 
from measurements of the azimuthal dependence of pion pairs. 

Note that we have used Zfykx, rather than /cy, as the transverse momentum variable, 
since this represents the transverse momentum of the pion pair relative to the jet. 

When Zakx gets large, the asymmetry decreases. This is expected, since the quark 
is then far off-shell, and the graphs approach their values with zero quark mass. When 
the quark mass is zero, there is exact helicity conservation along the quark lines, and 
hence there is no transverse-spin asymmetry. It can be easily checked from Eq. (4.10) and 
Eq. (3.5) that our calculated asymmetry is proportional to l/Zcy at large /cy. 

At small Zo, the calculated asymmetry appears to get smaller, with a broader distri- 
bution in ZfjkT- This may be a good prediction. But our model should not be reliable at 
small z. A reasonable prejudice is only that the spin dependence should decrease at small 
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z, since there one expects hadron production to be independent of the flavor and spin state 
of the initiating quark: The associated pomeron and gluon physics are not in our model. 



5. Conclusion 

Using the Linear Sigma Model to describe the fragmentation of a polarized quark has 
demonstrated the existence of a nonvanishing asymmetry from which polarization infor- 
mation can be obtained. Although the model should only be considered to give crude 
qualitative information, the large asymmetries we calculate do show that the spin correla- 
tion can be completely unsuppressed. 

This supports the value of doing experiments to measure the asymmetry experimen- 
tally. It can be used, for example, as a method of obtaining the transverse spin dependence 
of quarks in a transversely polarized proton.0 

It is clearly important to find better models for the fragmentation that include spin 
effects. Obviously, real fragmentation of a quark results in many more than two pions. 
This can easily dilute the effect of the physics we have modeled. Nevertheless, the number 
of pion pairs is fairly low if we restrict our attention to large z and fairly low invariant 
mass. 
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